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Abstract 

A covariant formula for conserved currents of energy, momentum and angular- 
momentum is derived from a general form of Noether’s theorem applied directly 
to the Einstein-Hilhert action of classical general relativity. Energy conservation 
in a closed hig-hang cosmology is discussed as a special case. Special care is 
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Is Energy Conserved in General Relativity? 

This may seem a strange question to be asking in 1997 considering that conservation 
of energy in general relativity was established by Einstein himself in 1916, the year 
after the theory was published [1]. Yet there remains a sour taste to Einstein’s pseudo¬ 
tensor formulation because of its lack of covariance. Modem attempts to improve on it 
can only really be said to be successful in specific cases such as asymptotically flat 
space-times and different forms of pseudo-tensor have only added to the ambiguity. 

To exemplify the problem it is helpful to ask again a question which is often posed by 
students. Where does the energy of the Cosmic Background Radiation go when it is 
red-shifted? Each photon carries a quantum of energy which slowly decreases as the 
universe expands. Ignoring interactions with matter which are irrelevant on the 
cosmological scale, the total number of photons remains constant. The total energy in 
the radiation must therefore decrease, where does it go? To avoid objections let us 
assume that the universe is closed so that the total number of photons is finite. 

The obvious answer is that the energy is transferred to the gravitational field which has 
an increasing (negative) total energy. But to justify this response we must be able to 
provide a general formula for energy which can be applied to this case. The pseudo¬ 
tensor is not sufficient because it is co-ordinate dependent and it is not possible to 
define a co-ordinate system over a whole spatial hypersurface of the closed universe if 
its topology is that of a three-sphere. How can we then be sure that we have not 
missed a leak of energy at the inevitable co-ordinate singularity where the formulation 
breaks down? 

It may well be possible to stop the gap with some mathematical slog, but a survey of 
most of the popular textbooks on general relativity does not turn up such an analysis. 
Those which use modem differential geometry shy away from the pseudo-tensor 
altogether and often leave general energy-momentum conservation outside the scope 
of the book. This is surprising because there is an elementary covariant way to treat 
energy and momentum conservation as I will show in this paper. Aside from the 
example I have given, the problem is important because of other cases such as energy- 
momentum transferred by gravitational waves and the contribution of gravitational self 
energy to mass which has been verified in tests of the strong equivalence principle by 
laser ranging of the moon. 

An interesting twist to the problem of total energy in a closed universe is that it comes 
out to be exactly zero rather than an arbitrary constant. This fact is often quoted when 
people ask where all the energy came from to set off the big bang, the answer: There is 
none! Again, it would be useful to have a tidy way of proving this. Sometimes the null 
total energy is taken as an indication that energy conservation in general relativity is an 
empty statement. This is not true because the energy is only zero dynamically, not 
kinematically. I.e. it is only zero when the gravitational field equations are satisfied. To 
see this we must provide the right formula for energy which can be applied to 
kinematics. 


2 



Energy Momentum in GR 


Energy and Noether’s Theorem 

Let us turn now to the gravitational field equations, 

R>‘'' = %nGT>‘'' 

By the Bianchi Identities, the divergence of the left hand side vanishes and implies, 

= 0 

If we were doing special relativity the divergence of the energy-momentum stress 
tensor would be adequate to establish energy-momentum conservation, but in general 
relativity it is not unless there is a killing vector field K , which generates an isometry 
of the metric. 


= K^T/ = 0 

When we have such a divergenceless vector current for energy flow, we can derive 
conservation because. 


1 d4gr 

' rs dx' 

and therefore the current can be integrated over a closed space-like hypersurface of 
constant time co-ordinate to give a conserved Energy, 

£ = f 

When we have no Killing field and we try this with the energy-momentum tensor we 
get. 



The second term, which can not be integrated, signifies a flow of energy between 
matter and the gravitational field. A term representing a contribution to the energy 
momentum tensor from the gravitational field is required. In deriving it Albert Einstein 
got some help from Emmy Noether who provided a general theorem relating 
conservation laws to symmetries of Eagrangians [2]. Noether’s theorem could be 
applied to the Hilbert-Einstein action for gravitation after an integration by parts 
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designed to remove second derivatives of the metric tensor from the action. This leads 
to Einstein’s non-covariant energy-momentum pseudo-tensor for which, 

A concise exposition of the derivation can be found in Dirac’s book on General 
Relativity [3]. After further manipulations which are described in the book of Landau 
and Lifshitz [4], a version of the pseudo-tensor which is symmetric in the two indices 
can be found. This symmetry is important if conservation of angular momentum is also 
to be demonstrated, but it makes the analysis even more complicated and less well 
motivated since Noether’s theorem is not used directly. Many other forms of the 
pseudo-tensor are possible such as the symmetric one given by Weinberg [5]. Recently, 
Bak, Cangemi and Jackiw derived a symmetric pseudo-tensor using Noether’s theorem 
and provided a useful comparison with other pseudo tensors [6]. Their result is still not 
covariant and is presented in a form which is divergenceless even kinematically. The 
full formula for the energy which shows the contribution to energy from different 
sources, and which is conserved dynamically but not kinetically, is often simplified 
using the field equations. 

With so many different versions of energy-momentum conservation which do not 
give exactly the same answer, one might ask what energy really is. The different 
formulations all differ by a quantity which is kinematically conserved (perhaps even 
before the simplification using the field equations) and which is asymptotically 
negligible in the weak field approximation to gravitational waves. Even so, we might 
have expected a unique concept of energy. 

The correct answer to the question “What is energy?” is that the correct energy is that 
which is given by Noether’s theorem. This is still not unique since the Lagrangian is 
only unique up to a term which is a perfect divergence. This feature is what is used to 
get rid of the terms with second derivatives in the metric tensor before deriving the 
energy equation, but it also destroys covariance. Even with covariance other forms of 
the action are known but at least it might be worth while having the correct expression 
for energy for the traditional Einstein-Hilbert action. 


16 ;rGJ 

This can be done using a form of Noether’s theorem which works with second 
derivatives in the Lagrangian. 


(-4sR + LM)d^x 
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Consider a general action of the form, 

S=jL(f 

(p' are all the components of the field variables including the metric tensor. Applying a 
variation principle, the Euler-Lagrange equations are. 


f 


dL . 


dL 




f 


+ (9 






dL 




(90' 


= 0 


J 


A continuous symmetry of the action which might include co-ordinate changes as well 
as field transformations would have a general infinitesimal form, 

00“ = £0“ 

5x^ = ee 

and invariance implies that, 

SL = -^ee‘‘ +^^ee\^,=d (ek^L) 

Combining results we finally derive the conserved Noether current. 


1^ = 


dL 

dr. 
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= 0 
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The Covariant Formulation 

The next step is to apply this form of Noether’s theorem directly to the Einstein- 
Hilbert action. The action is invariant under diffeomorphisms (co-ordinate 
transformations). Infinitesimal diffeomorphisms are generated by 

= rat" 

Sgaf, = e^8al,=e(K,f,+’^f,J 
y/‘ =eff(k>‘) 


Ip' are all the field components other than the metric tensor. The answer will not take 
the form of an energy-momentum tensor, it will be a contravariant vector current 
which depends on the contravariant transport vector field k and by Noether’s theorem 

a 

it will be conserved for any choice of k . This is a natural consequence of the fact that 
in general relativity the symmetry is diffeomorphism invariance rather than global 
Lorentz invariance. 


The derivation is quite long but will nevertheless by given in some detail for the 
gravitational part. Our starting point is the Lagrangian density, 


L = 


-1 

16;rG 




It is assumed that the matter Lagrangian contains no higher than first covariant 
derivatives in field variables, and transforms as a scalar-density under diffeomorphisms 
In that case the matter field contribution to Noether current is given by, 

=J^G + J^M 
J^M = 


The Riemann curvature tensor is given by. 


= R 


— 


R 


gvpa 
( 1 ) 


jivpo 


+ R 


( 2 ) 


pvpa 


pvpa 


1 ^§pG,vp Sva ,pp §pp,va §vp ,fj.a ) 


,lip o fip ,va o vp ,pa ^ 
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R = Rpvpag 8 
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The gravitational contribution to the current will be split up accordingly. The factor 
passes outside the partial derivatives but an extra term must be included when it 
passes through the co-ordinate derivatives. 


J>‘G=yk;{k^R-f) 






^Sap ~ ^a-p ^p-a 


(dy +r VCt) 


dR 


( 1 ) 


A 






^^aP ^^Sap\v 


dg 


ap,jxv 


and we can calculate, 


3p(i) 

_ gl^iP gCC)V 

O O O O 

SaP,jxv 

If we knew that the answer must be covariant, then we could use the fact that, in a 
local inertial frame, the first derivatives of the metric tensor vanish at a single point. 
This makes all the difficult terms in the current vanish and the calculation is easy to 
complete. 


However, covariance has not yet been proven, so the terms should be worked out at 
length to be sure of the result. 
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-/"c = i6fc (<:"«-2r„'‘+r;'‘„ +r “„) 


This is the final covariant expression for the contribution to energy-momentum current 
from the gravitational field. The dependency on the contravariant transport vector k is 
left explicit. Any attempt to remove it to derive a stress tensor will destroy covariance 
and make the result co-ordinate dependent. The generality of the expression would be 
lost. The term with the curvature scalar depends on second derivatives of the metric 
field in addition to first derivatives. This is unusual for an expression of energy or 
momentum. The second derivatives could be removed by adding a term which is 
kinematically divergence free but only at the cost of covariance. 

The choice of transport vector field is used to distinguish between currents of energy, 
momentum, angular-momentum, etc. In the case where it is everywhere time-like the 
current is interpreted as flow of energy. Energy and momentum are therefore not 
unique but this is no surprise or obstacle. For a given co-ordinate system the energy 
and momentum which are conjugate to time and space co-ordinates can be found but 
they do not form a tensor. 

Using the identity. 





and including the contribution from the matter fields we get 
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J" = ukG(k'‘R- +k^“a- 2k“R /) + k“T/ 

If we now apply the gravitational field equations the expression can be simplified 
dramatically. 


~ i67iG ^ 

This can also be written in a form without explicitly covariant derivatives, 

The divergence of this expression is identically zero but it is only correct dynamically. 
The full expressions we derived before simplifying with the field equations are 
kinematic formulae for the energy-momentum currents and their divergence is not 
identically zero, but it is dynamically zero. They show directly the contributions to 
energy and momentum from separate parts of the Lagrangian. 

These final dynamic expressions are the Komar superpotential [7]. When applied to an 
asymptotically flat space-time with appropriate use of boundary conditions this leads 
on to the Geroch-Winicour linkage formulation [8], the Bondi Energy and the ADM 
Energy [9]. 

Energy in a Closed Cosmology 

We are now in a position to consider the case of energy in a closed cosmology. The 
covariant currents we have constructed are dynamically divergenceless. When 
integrated over a closed space-like hypersurface to give the total energy or momentum 
of the universe, the result must be zero because there is no boundary. As we noted in 
the introduction, this is hard to prove rigorously with a pseudo-tensor formulation, but 
with covariant currents it is a simple consequence of Gauss’ theorem. In a closed 
universe with 3-sphere spatial topology there is no clear distinction globally between 
momentum and angular-momentum. They are consequences of rotations of the spatial 
3-sphere about different points. The Minkowski model of energy-momentum 4-vectors 
and angular-momentum tensor is therefore not relevant to a closed spherical spatial 
topology. 

a 

Eor energy we must be careful about our choice of the transport vector field k . If its 
covariant form k« had components which were constant in some co-ordinate system 
then the current would be zero everywhere. That would not necessarily be incorrect 
but it should normally be chosen so that it increments the time co-ordinate uniformly. 

It follows that the energy current depends critically on the choice of time co-ordinate 
as is to be expected. 

In a homogeneous Eriedmann-Robertson-Walker cosmology there is a natural choice 
of time co-ordinate, that of the co-moving co-ordinate system. In a more realistic 
locally inhomogeneous cosmology the matter is not so simple. We could define 
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cosmological time at each space-time event to be the maximum proper-time integrated 
along all time-like paths from the big-bang singularity to the event. This is well defined 
provided the cosmology has no closed time-like curves and can be sequentially foliated 
into space-like surfaces. It is easy to show that the surfaces of equal cosmological time 
are continuous space-like surfaces, except where they meet singularities which have a 
finite age. Even without singularities such surfaces are not smooth because co-ordinate 
systems which incorporate cosmological time are synchronous reference frames. It is 
known that they tend to develop co-ordinate singularities in realistic space-times (see 
[4]). In general, the energy currents can be defined relative to a smooth space-like 
foliation of space-time and a time-like transport vector field It is beyond the scope of 
this article to establish when they exist. 

Black-hole singularities do not actually cause a breakdown of energy-momentum 
conservation because it is possible to pass a space-like surface through a black hole 
which avoids the singularities. This can be done even at late times. The surface would 
have to bend back in time at nearly the speed of light but can narrowly avoid 
becoming time-like. This can be seen in a Carter-Penrose diagram of a black-hole. 
Fortunately the total value of energy momentum and angular momentum for an 
isolated object depends dynamically only on the choice of transport vector at the 
boundary of the surface and not on the surface or choice of transport vector internally. 
This is not true kinematically. 

Because of the difficulty of choosing a time co-ordinate which can be applied 
universally and which avoids all singularities, the interpretation of total energy in a 
closed universe is not perfectly clear. This is not too surprising or worrying since 
energy and momentum are usually considered to be relative, and the whole universe 
cannot be considered relative to anything else. Nevertheless, there are covariantly 
conserved energy and momentum currents whose total integrates to zero. This is not a 
trivial result since it is not true kinematically if we use the correct kinematic definition 
of the currents which is derived directly from Noether’s theorem without using the 
field equations to simplify. Remember that it is only this kinematic form which shows 
explicitly the contribution to the currents from individual fields. 

Conclusions 

I have shown how to derive conserved energy and momentum currents by applying 
Noether’s theorem directly to the Hilbert-Einstein action of general relativity. I have 
not entered into questions of lacalisability of the gravitational energy and have not 
considered how to choose the transport vector or boundary conditions to obtain 
physically useful answers in particular cases except for a closed cosmology. 

Note Added 

I am informed that the Komar superpotential [7] may have been derived from a 
variational principle in this way by Trautman in a conference paper in 1961. A more 
detailed derivation was given by Kijowski in 1978 [11]. It is my pleasure to thank Piotr 
Chrusciel, Joseph Katz, Jerzy Kijowski, Greg Oganessyan, and Kumar Shwetketu 
Virbhadra for these and other useful references and comments. 
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